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We develop a nonlinear theory of soliton-induced waveguides that describes a finite-amplitude 
probe beam guided by a spatial dark soliton, in a saturable nonlinear medium. We suggest an 
effective way of controlling the interaction of these soliton-induced waveguides and also show that, 
in sharp contrast with scalar dark solitons, the dark-soliton waveguides can attract each other and 
even form stationary bound states. 


It was suggested in the early days of nonlinear optics, 
that the waveguiding properties of a self-trapped laser 
beam could be used to guide atoms or ions along the 
beam’s axis Q. A similar physical mechanism is respon¬ 
sible for the trapping of a weak beam by a spatial (bright 
or dark) optical soliton. This phenomenon is usually re¬ 
ferred to as the soliton-induced waveguiding of a probe 
beam. It has been widely discussed theoretically and has 
also been observed experimentally 

As has already been established, the major advantage 
of employing the waveguiding properties of dark, rather 
than bright, spatial solitons is their greater stability and 
steerability in a Kerr medium |^. Moreover, in the case 
of saturable nonlinearity, dark-soliton waveguides remain 
single-moded , unlike the waveguides created by bright 
solitons 1^. 

The theory of soliton-induced waveguides developed so 
far is based on the assumption that a probe beam guided 
by a spatial soliton is weak (see, e.g., Refs. |^,^,^). How¬ 
ever, when the intensity of the probe beam is increased, 
the guided mode requires a broader waveguide and the 
beam interacts strongly with the waveguide, changing it 
in a self-consistent manner. Such a structure can not be 
described by a linear waveguide theory because it rep¬ 
resents a bound state of the coupled bright and dark 
components. In this Letter we develop a nonlinear the¬ 
ory of the soliton-induced waveguides, treating the case 
of two incoherently interacting beams in a photorefrac- 
tive medium as a model already realized experimentally 
1^. We also analyze the interaction of two neighboring 
waveguides and show that attractive forces acting be¬ 
tween the guided beams can result in suppression of the 
repulsion between dark solitons. 

Following the original paper by Christodoulides et al. 
, we consider two incoherently interacting linearly po¬ 
larized beams with the normalized envelopes U and W 
in a photorefractive medium, described by the equations, 
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Here p = loa/Id, /? = {koXo)‘^n‘^r33Eo/2, where loo 
stands for the total power density away from the beam, 
Id is the so-called dark irradiance, fco is the propaga¬ 
tion constant, xq is the spatial width of the beam, and 


^e^ 33^'0 is a correction to the refractive index that it due 
to the external field applied to a crystal along the x-axis 
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Localized solutions of Eqs. (|1|) can be sought in the 
form U = and W = 

and they form a continuous two-parameter family. Sys¬ 
tem can be further simplified by measuring the spa¬ 
tial coordinates z and x in units of •\/s//3(l -|- p) and 
{s/[/3(l -I- respectively, where s = 1 — gi. Then, 

u and w satisfy the normalized equations: 

du 1 d'^u (ImP -I- |wP)rt 
dz 2 dx^ 1 + s(|m|2 + |u;|2) ^ “ 

.dw 1 d'^w (|up-I-Iwpjrc 
* 9z 2 dx^ 1 -I- s(|u|2 -I- |w|2) ^ 

where A = {1 — q 2 )/s < 1 is the dimensionless soliton 
parameter, and s characterizes nonlmearity saturation, 
(s < 1). In the limit s —> 0, Eqs. (||) reduce to the 
Manakov model which possesses exact bright and dark 
multisoliton solutions [pTHl^. 

Stationary solutions are described by Eqs. (||) with 
z derivatives omitted. We assume that the component 
u has nonvanishing asymptotics, whereas w is spatially 
localized. Then the simplest solution of this kind de¬ 
scribes a dark soliton of the field u, provided that w = 0. 
This solution can be presented in quadratures and then 
found numerically |^. For the component w, the dark 
soliton creates an effective waveguide. To analyze what 
kind of guided mode can be supported by this waveguide 
we linearize the equation for w and study its localized 
solutions. That is, we consider a one-component dark 
soliton u = Us(x) and add a small perturbation written 
as w = 0(e) and u = Us{x) + O(e^). It is easy to verify 
from Eqs. (^) that 0(e) in the expansion for u is zero. 
Substituting this expansion into Eqs. (|^), we obtain two 
decoupled equations, of which the stationary equation for 
w is 
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w + Xw = 0. 


( 3 ) 


Depending on the properties of the effective waveguide 
V{x) = |Ms(a:)p(l-I-s|us(x)p)“^, eigenvalue problem (||) 
possesses a number of solutions that decay as a; —> ±oo, 
the so-called guided or bound modes. In the limit s ^ 0, 
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the symmetric solution ryi(a;) = sechcc exists for any 
A > Ao = 1/2. In general, the cutoff of a bound mode, 
Ao(s), depends on the saturation parameter s. Further¬ 
more, if Us{x) is a dark soliton of a saturable medium, 
there exists only one bound mode of eigenvalue problem 

t for any value of s < 1, as was also found in Ref. 

, Therefore the dark-soliton waveguide in a saturable 
medium is always single moded . 



FIG. 1. Bifurcation of a one-component dark soliton Ub(x) 
for three values of s, shown as the dependence of the compli¬ 
mentary power Pc on the propagation constant A. Examples 
of two-component localized solutions are shown for s = 0.2. 



FIG. 2. Evolution of the dark-soliton waveguide without 
the guided mode for s = 0.5. (a) Soliton profiles at z = 0 
(thick line) and z = 100 (thin line) for A = 0.95. (b) Ampli¬ 
tude of the additional gray-soliton pair g versus A measured 
at 2 = 100. 
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FIG. 3. Examples of two-soliton bound states in model (|^: 
(a) s = 0.4 and A = 0.55, (b) s = 0.5 and A = 0.55, (c) s = 0.5 
and A = 0.65. 

The localized solution that describes the fundamental 
(symmetric) mode of a dark-soliton waveguide appears 
as a bifurcation of the one-component dark soliton, Us, 
in which a two-component solution (u, w) emerges from 
the one-component solution (us,0) at certain A = Ao(s). 
Such localized solutions form a continuous family in A 
for a fixed s. Figure 1 presents examples of such a bifur¬ 
cation together with the corresponding two-component 


soliton states, found numerically by the relaxation tech¬ 
nique, for different values of s. Shown is the dependence 
of the complimentary power, 

/ -l-oo 

dx{\u{x)\^ + \w{x)f - (1 - s)“^} , 
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on the soliton parameter, A. 

From the physical point of view, the parameter A char¬ 
acterizes the intensity of a bright component guided by a 
dark-soliton waveguide. As A grows, the intensity of the 
bright component grows as well, remaining less than that 
of the dark-soliton background. As it approaches a max¬ 
imum, the dark soliton becomes broader. This means 
that, for A much larger than Ao(s), the soliton waveguide 
deforms: It becomes wider to guide a beam of higher 
intensity. This effect is definitely beyond the linear ap¬ 
proximation provided by the model [Eq. (^]. It can, 
however, be described by a nonlinear theory constructed 
by numerical solution of coupled equations (||). The cor¬ 
responding solutions, presented in Fig. 1 for A > Ao(s), 
generalize the localized solutions of the linear eigenvalue 
problem [Eq. (^] valid only near the bifurcation points 
A where A ~ Aq. 

To verify the concept of a deformed waveguide, we per¬ 
formed simulations of the waveguide propagation in the 
case when the bright component is removed. As expected 
from the theory of dark solitons 0 and previously con¬ 
firmed in 0, a wider dip in the input beam generates at 
least one additional pair of gray solitons, as is shown in 
Fig. 2(a). The amplitude of the additional gray-soliton 
pair g is shown in Fig. 2(b) as function of the parameter 
A. 

Similar to the case of the exactly integrable Manakov 
system ||^,|IJ] , the nonintegrable model (||) supports sta¬ 
tionary bound states of two (or more) solitons. We have 
found such stationary structures numerically. Examples 
of the stationary solutions, corresponding to different val¬ 
ues of A and s, are shown in Figs. 3 (a,b,c). 

The stationary bound states are expected when the to¬ 
tal force acting between the neighboring two-component 
solitons vanishes. Two closely separated dark solitons al¬ 
ways repel each other, and this property does not depend 
on integrability |^]. In the light of this knowledge, the 
existence of stationary solutions in the form of two bound 
dark-soliton waveguides is rather surprising. Therefore, 
we are led to wonder whether the introduction of a bright 
component causes an attractive force between solitons 
that could nullify the repulsion of the dark components. 
Note that the stationary bound states of the model [Eqs. 
(^)] are formed by bright components that are out of 
phase with each other (Figs. 3). Expecting them to at¬ 
tract each other would be inconsistent with the theory 
of bright solitons that predicts attraction only between 
in-phase solitons. However, in our case the bright compo¬ 
nent by itself is not a soliton because, after the dark com¬ 
ponent is removed, it diffracts in the defocusing medium. 
Localization of the bright components is due to trapping 
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in an effective waveguide created by dark components. 
As a result, interaction of these guided modes is different 
to that expected for bright solitons. From the physi¬ 
cal point of view, even slight overlap between two out 
of phase localized beams in a defocusing medium leads 
to an increase of the refractive index in the overlapping 
region, therefore increasing the attraction of the beams. 



FIG. 4. Interaction between two dark-soliton waveguides 
guiding bright components (a) in phase or (b) out of phase. 

Analyzing interaction between the soliton waveguides, 
we look for a solution to describe two weakly overlap¬ 
ping two-component solitons in the form: u = tanh(a; — 
xq) tanh(a:: -|- xq) and w = a sech(a: — xq) -I- a e*‘^sech(a: -|- 
cco). Then we employ the Lagrangian formalism for dark 
solitons jl^] to derive effective equations for the soliton 
parameters, including the soliton separation 2xq and the 
relative phase between the bright components (p. Omit¬ 
ting lengthy details, we present the effective energy of 
the soliton interaction in the form (a:o 3> 1), 


(xo) = + 4a^ a:oe“2““ cos(</)). (4) 


The first term in Eq. (|^) describes the interaction of 
two scalar dark solitons; see Eq. (34) in Ref. [^. The 
second term appears as the result of the interaction be¬ 
tween the bright components of the amplitude a. As 
follows from Eq. (^), the repulsive force between dark 
components is weaker than the force acting between the 
bright components. Moreover, introducing bright com¬ 
ponents into the closely separated dark solitons can lead 
to a nontrivial effect when the repulsion of the dark com¬ 
ponents is compensated for by an attractive force acting 
between the bright components. This effect requires that 
the bright components be out of phase (^ = tt) , in agree¬ 
ment with the stationary solutions found numerically; 
see Fig. 3. Figure 4 shows the results of the propa¬ 
gation of two neighboring dark-soliton waveguides when 
the bound modes they guide are in phase [Fig. 4(a)] or 
out of phase [Fig 4(b)]. In the former case the inter¬ 
action between bright components is repulsive, and the 
solitons repel each other even more strongly than two 
scalar dark solitons. In the latter case, the interaction 
is attractive and it forces the dark solitons to collide. 
These results confirm that the interaction of the dark- 
soliton waveguides is phase sensitive. We believe that the 


demonstrated behaviour can be observed experimentally 
for the incoherent interaction of photorefractive solitons. 

In conclusion, we have developed a nonlinear theory of 
the soliton-induced waveguides created by dark solitons 
in a saturable optical medium. We have found families 
of the soliton states describing a finite-amplitude guided 
mode trapped by a dark soliton. We have described a 
novel type of interaction of the dark-soliton waveguides 
in which the repulsion of the neighboring dark solitons 
is suppressed by an attractive force acting between the 
guided modes, leading to the existence of stationary two- 
soliton bound states. 
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